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Abstract. We present the first detection of the envelope which surrounds the post-AGB binary source HR 4049. VINCI-VLTI 
A"-band interferometric observations of this source imply the existence of a large structure with a Gaussian angular FWHM 
22.4 ± 1.4 mas or uniform disk diameter of 34.9 ± 1.9 mas. With the Hipparcos parallax of 1.50 ± 0.64 mas these values 
correspond to a physical size of 14.9_;V AU and 23.3*1)? AU, respectively. Our measurements, covering an azimuth range of 
~ 60° for the sky-projected baseline, provide information on the geometry of the emitting region and show that there is only 
a slight variation of the measured angular values along the different directions sampled. Thus, our results are consistent with a 
spherical geometry of the envelope. However, we cannot completely rule out the existence of an asymmetric envelope (like the 
circumbinary disk envisaged by some recent models) because of the limited spatial frequency and azimuth range covered by 
the observations. 
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1. Introduction 

HR 4049 (HD 89353, AG Ant) is considered the prototype 
of a peculiar group of sources in the post-Asymptotic Giant 
Branch (AGB) stage of evolution that are believed to be close 
binary systems surrounded by a massive and stable circumbi- 
nary envelope as a consequence of mass loss processes dur- 
ing the AGB phase. These sources typically show evidence of 
extremely metal-depleted photospheres, possibly the result of 
dust formation in the circumstellar environment followed by 
re-accretion of the depleted gas onto the stars. In the case of 
HR 4049 the presence of large amounts of circumstellar dust 
is indicated both by the UV deficiency and strong IR excess 
( ILamers et al. 198fjl . At the moment, however, the size and 
spatial distribution of this material are not well known. 

Several models of HR 4049 have been proposed, mainly 
based on photometric and polarimetric observations, with dif- 
ferent predictions on the geometry and extent of the dust en- 
velope, depending also on the physical properties of the dust 
grains. Waters et al. (1989) suggested, for instance, that the 
observed IR excess could be modelled by the emission from 
an optically thin cloud with a radius ranging from 5 to about 
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300 AU or alternatively by a blackbody-emitting structure with 
a size of the order of 2 AU. More recent work on HR 4049 
tended to favour the flattened or asymmetric disk geometry 
for the distribution of the dust around the source. Photometric 
variations in phase with variations of radial velocity indeed 
suggest that the circumstellar matter is distributed in a disk 
surrounding the binary, which would be observed from Earth 
at a high inclination with respect to the plane of the sky 
(IWaelkens et al. 199U . Analysis of the spectral energy distri- 
bution (SED) of HR 4049 has shown that the infrared excess 
between 1 and 850 /mi can be fitted almost perfectly with 
a single blackbody function at T ~ 1150 K. On this basis, 
it has been suggested that the radiation we observe is com- 
ing from the inner rim of a circumbinary disk heated by the 
central sources ( Domin ik et al. 20 03 1. Moreover, polarization 
measurements indicate the presence of Rayleigh scattering by 
two different populations of grains, suggesting that the distri- 
bution of the circumstellar material might be asymmetric with 
the presence of a disk and a bipolar structure ( J oshi et al. 19871 
Johnson & Jones 1991 Johnson et al. 1999 1. 

In this paper, we present the analysis of VLTI- VINCI ob- 
servations intended to shed light on this interesting problem 
by measuring the size and shape of the circumstellar enve- 
lope of HR 4049 directly by interferometric means. The ob- 
servations clearly resolve for the first time a relatively large 
structure around the source with a dimension and shape that is 
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Table 1. Observational log with nights, baselines and sky- 
projected baseline azimuth for VINCI observations of HR 
4049. 



Night 


VLTI Baseline 


Proj. Baseline (m) 


Azimuth 


2002-12-16 


B3-C3 


8m 


7.4-7.5 


53°-62° 


2003-01-13 


B3-C3 


8m 


7.8-7.9 


75°-76° 


2003-01-15 


B3-C3 


8m 


7.8-8.0 


64°-76° 


2003-01-16 


B3-C3 


8m 


7.9 


74°-75° 


2003-01-20 


B3-C3 


8m 


7.5-7.8 


77°-80° 


2003-02-25 


B3-D1 


24 m 


11.3-12.3 


105°-108° 



quite plausibly attributable to the putative circumbinary enve- 
lope surrounding the binary. 

2. Observations and data reduction 

The observations (see Tab^for a summary) have been carried 
out with the VLT Interferometer (e.g. Glindema nn et al. 2 003 1 
between December 2002 and February 2003 using the near- 
infrared /T-band VLTI commissioning instrument VINCI (e.g. 
IKervella et al. 200*31 . which recombines the light coming from 
two 0.35 m aperture siderostats (IDerie et al. 20 03 1. The data 
were taken in the framework of the VLTI shared risk sci- 
ence program and are publicly available from the ESO Science 
Archive. The two interferometric baselines B3-C3 and B3-D1 
were used, which have a ground length of 8 and 24 meters, re- 
spectively, and the same ground position angle (East of North) 
of71°. 

With this configuration, several series of interferograms 
(consisting of 100, 200 and 500 scans) were obtained on HR 
4049, with a fringe frequency of 60 Hz. This rather slow fre- 
quency had to be chosen because of the relatively low K- 
magnitude of the object (K = 3.21 + 0.29. ICutri et al. 20031 
The details about the calibrator stars used for the observations 
are shown in Tab|5] The fringe frequency for several calibra- 
tors was 290 Hz, different from the one used on the science 
target. This may introduce a calibration error which has been 
taken into account. By comparing the results obtained using 
only calibrators sampled with the same fringe frequency as the 
target (whenever possible) to those derived using all available 
calibrators, we have estimated the effect of this error on the 
calibrated visibilities to be less than 3%. 

Squared coherence factors were calculated processing the 
raw interferograms with the VINCI data reduction software 
(version 3.0) by Kervella et al. (2004). Subsequent calibration 
of the data and calculation of broad-band synthetic visibility 
values were performed as in Wittkowski et al. (2004). 

A total of 17 values for the squared visibility were finally 
obtained from data reduction at an effective wavelength of 
/left = 2.21 fim and in an azimuth range for the sky -projected 
baseline between 53° and 108°(North through East). 

We use two simple two-component models for the object 
intensity distribution, both composed of a single central star 
(the primary of the binary system) sourrounded by an envelope. 
For the central star, we use the radius R* = A-IRq derived by 
Bakker et al. (1998) and adopt an almost unresolved uniform 
disk (UD) brightness distribution with angular diameter 0.7 



Table 2. Calibrator stars. The uniform disk diameter #ud and 
the corrisponding error cr g are reported together with spectral 
type, K magnitude and effective temperature. 



Star 


i3U. ± vpc 


Kmag 


"ID 


CT ft 
















31 Ori" 


K5 III 


0.90 


3.56 


0.057 


4046 


51 Ori" 


Kl III 


2.10 


1.87 


0.021 


4508 


a Hya'' 


K3II-III 


-1.21 


9.11 


0.911 


4230 


6Lsp a 


K0III 


1.37 


2.56 


0.041 


4656 


eLep° 


K5 IIIv 


-0.20 


5.91 


0.064 


4046 


i Hya" 


K2.5 III 


0.91 


3.41 


0.048 


4318 


v 2 CMa" 


K1.5 III/IV 


1.62 


2.38 


0.026 


4497 


yuHya 6 


K4III 


0.40 


4.55 


0.455 


4090 


v337 Car" 


K2.5 II 


0.03 


5.09 


0.058 


4300 



" from Borde et al., 2002, A catalogue of calibrator stars for LBSI. 
6 from calibration by Dyck et al., 1996 (Radii and Effective 
Temperatures for K and M Giants and Supergiants), using K magni- 
tudes from Gezari, 1999 (Catalog of IR observations). 

mas (on the basis of the Hipparcos-determined distance). Since 
the flux coming from the source is mainly due to the primary 
(e.g. Bak ker etal. 1 998) we do not take into account the much 
fainter secondary which has been detected by radial velocity 
measurements ( Waelk ens et al. 199U IBakker et al. 1 998). As 
for the envelope, we describe it with an UD and a Gaussian 
distribution. We assume a ratio F = F t /F tot ~ 15% for the stel- 
lar flux with respect to the total observed flux in the ^T-band 
JDominik et al. 20031 . Keeping the adopted parameters of the 
(unresolved) star constant and treating the UD diameter or the 
Gaussian FWHM of the envelope as the only free parameter, 
we compare the observed visibility with our simple model. 

Although these simple two-component models certainly 
represent a first approximation, this procedure can provide us 
basic information on the shape and dimensions of the emitting 
region in HR 4049 without having to resort to more complex 
and uncertain assumptions. We should bear in mind that inter- 
ferometer observations carried out at different azimuth values 
for the projected baseline allow us, in principle, to detect asym- 
metries of the intensity distribution. 

For this reason, we grouped the data points in four differ- 
ent 10° wide azimuth bins (53°-62°, 63°-72°, 73°-82° and 
103°-112°). Then, we performed a best-fit of our models to 
the squared visibility values for each bin to obtain the corre- 
sponding values of the Gaussian FWHM and UD diameter of 
the envelope (see Fig.^and Tab[3}- 

We also investigated how the obtained results vary depend- 
ing on the adopted flux ratio F. We found that, when changing 
F by 5%, the related variation for the envelope angular size is 
small: about 5% and 4% for the Gaussian and UD distribution, 
respectively. 

Finally, although the decrease of the observed visibility be- 
tween spatial frequencies ~ 16 and ~ 26 cycles/arcsec clearly 
shows that we are resolving a compact envelope component, 
a contribution to the total visibility by an additional extended 
scattering component is quite possible. In this case, assuming 
that 10% of the total flux is due to this incoherent compo- 
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Table 3. VLTI- VINCI measurements of the size of the enve- 
lope surrounding HR 4049. Depending on the adopted source 
model (UD + Gauss or UD + UD, see the text), the FWHM of 
a Gaussian distribution or an equivalent uniform disk diameter 
6 ™ are reported. Formal errors are shown; further error (~3%) 
is introduced by the calibration process (see the text). The mea- 
surements refer to the four azimuth bins mentioned in the paper 
and to a value F—15% of stellar to total flux ratio. The last row 
shows the results when considering all data together. 



Azimuth 


UD + Gauss 


UD+UD 


Bins 


FWHM ±cr (mas) 


0^±a-(mas) 


53°-62° 


22.7 ±0.1 


36.9 ±0.1 


63°-72° 


23.3 ± 0.4 


37.7 ± 0.6 


73°-82° 


23.8 ± 0.6 


38.4 ±0.8 


103°-112° 


21.7 ± 0.2 


33.4 ± 0.4 


all azimuths 


22.4 ± 0.5 


34.9 ± 0.7 



nent, we found that the fitted sizes of the envelope (UD and 
Gaussian) decrease by roughly 10%. 

The measurements show that we have detected an emitting 
structure whose azimuthally averaged angular size is 22.4 + 
1.4 mas (Gaussian FWHM) or 34.9 + 1.9 mas (UD), with er- 
rors including also the previously mentioned uncertainties due 
to calibration and to the adopted flux ratio. There is only a 
slight variation of the size depending on baseline orientation, 
most of it when passing from the first three adjacent bins to 
the fourth one; the values range from 21.7 to 23.8 mas for 
the Gaussian and from 33.4 to 38.4 mas for the UD distribu- 
tion. This indicates that our observations are consistent with a 
roughly spherical symmetric envelope; however, having mea- 
surements for only two different baselines, and therefore in a 
limited range of spatial frequencies and position angles, does 
not allow us to definitely rule out at this point the possibility 
that the envelope might be partly asymmetric and hence to un- 
ambiguously discriminate between a disk structure or a more 
spherical geometry. A more flattened envelope could explain 
the small variability of the object in the /f-band (about 0.1 mag 
in time scales of the order of 250 days, Waelke ns et al. 19 91 1; 
this would be due to different extiction along varying lines 
of sight to the primary as the components orbit each other 
(Waelken s et al. 199U . Because of our limited data sample we 
can only give here an estimate of the envelope size and we can- 
not exclude the existence of more complex structures consist- 
ing of several clouds or regions with different optical depth. 

With the Hipparcos parallax of 1.50 + 0.64 mas 
jPerryman et al. 199 7 ) the azimuthally averaged angular values 
obtained for the Gaussian FWHM and the UD diameter of the 
envelope (Tab|3} correspond to linear dimensions of I4.9+44 1 
AU and 23.3^ 7 3 AU, respectively. 

3. Discussion 

The geometry and extension shown by our results are in agree- 
ment with those proposed by Waters et al. (1987), who mod- 
elled the emission from HR 4049 as arising from a simple 
spherical optically thin envelope and found that such structure 
would have a radius ranging from 5 AU up to about 300 AU. 



Joshi et al. (1987) suggested that the observed intrinsic po- 
larization of light from HR 4049 is introduced by scattering due 
to the circumstellar matter and that the low degree of polariza- 
tion measured probably indicates that the envelope is almost 
spherically symmetric, which would be in agreement with our 
results. They modelled the strong dependence of the position 
angle with wavelength assuming the presence of two spheroidal 
dust shells surrounding the source, having different symme- 
try axis and composed by two different populations of large 
(a ~ ljum) and smaller grains (a < 0.05//m). 

Johnson and Jones (1991) modelled the polarization as aris- 
ing from scattering by dust in an envelope having the shape of 
a prolate ellipsoid; they found a position angle of about 8 de- 
grees, remarking how this angle is approximately perpendicu- 
lar to the direction of an extended (about 3 pc) emission at 100 
jum reported by Waters et al. (1989), presumably due to cold 
matter ejected from the source during the previous evolution- 
ary stage. 

It is also of interest to compare our interferometric results 
with the model proposed by Dominik et al. (2003) in which 
the emission we observe would actually come from the inner 
surface of a geometrically (Height /Radius ~ 1/3) and opti- 
cally thick circumbinary disk surrounding the central sources 
and seen at an inclination of ~ 60° with respect to the plane 
of the sky. Using this model, the authors infer a radius of or- 
der R ~ 10 AU for the distance of the inner surface of the disk 
from the star, which is roughly consistent with the size of the 
structure detected by VINCI. Because of the high inclination, 
however, we would expect to observe a somewhat larger vari- 
ation in the measured size of the structure for different base- 
line orientations. Assuming that we are observing in an azimuth 
range around the major axis of the projected disk (this would 
reproduce the trend we observe in our measurements, with the 
smallest values for the first and last bin), we should in fact see 
a diameter variation of at least 4 mas, which is around the 3cr 
level in our measurements. 

Finally, Johnson et al. (1999) proposed that the measured 
optical polarization indicates Rayleigh scattering arising in an 
optically and geometrically thin circumbinary disk. This, to- 
gether with varying extinction with orbital phase, has led the 
authors to formulate the hypothesis that the orientation of the 
disk is tied to the orbit of the binary and so the projected face 
of the disc would change as the central sources orbit each other. 
Moreover, the UV polarization shows a typical 90° change 
in position angle, usually related to the presence of a bipolar 
structure (i.e. two structures orientated perpendicularly to each 
other), which is often seen in planetary nebulae. In their mea- 
surements, the polarized light shows two preferential position 
angles (averaged over wavelength) of ~ 17° and ~ 28°. It is 
interesting to notice that, since the observed PA of scattered 
light is supposed to be perpendicular to the plane of the cir- 
cumstellar material in the disk, the orientation angle of the long 
axis of the disk, as viewed from Earth, would lie in the range 
1 10°-120°. Thus, the VINCI observations presented here (with 
a ~ 50°-110°) would sample directions ranging from approx- 
imately that of the (projected) disk's long axis (a ~ 110°) to 
one almost perpendicular. However, our results show that there 
is no clearly evident difference in the source dimensions when 
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passing from one azimuth bin to the others, as one might expect 
assuming the disk configuration reported above. 

Although HR 4049 is generally believed now to be sur- 
rounded by a circumbinary disk on the basis of the different 
observational evidences mentioned above, the existence of a 
more spherical envelope should not be excluded out of hand, 
as the recent example of the Herbig AeBe stars clearly shows 
(e.g. lMiUan-Gabet et al. 200TV 

Moreover, it is also likely that the envelope around this 
source might be multi-structured, as the polarization data sug- 
gest (bipolar structure, different populations of grains in sepa- 
rated shells). 

4. Conclusions 

The reduction and analysis of VINCI data on HR 4049 have led 
to the detection of a structure of angular dimensions 34.9 mas 
(uniform disk diameter) or 22.4 mas (FWHM of a Gaussian 
distribution), probably the circumstellar envelope or disk sur- 
rounding this binary system. 

Assuming the Hipparcos-determined distance of 667 pc for 
the source, this angular values corresponds to a size of about 
23.3 AU and 14.9 AU respectively. 

The measured size shows only a slight variation along 
different directions spanning an azimuth range of ~ 60°, 
consistent with models of spherically symmetric or ellip- 
soidal envelopes (e.g. IWaters et al. 19891 IJoshietal. 19871 
IJohnson & Jones 199U . 

We cannot, however, completely rule out the possibility of 
the presence of asymmetric structures, such as the circumbi- 
nary disk predicted by some recent models. 

Our knowledge of the geometry of HR 4049 system can 
be refined with further interferometric observations. Sampling 
different azimuths and spatial frequencies and also carrying out 
measurements at different wavelengths would definitely help 
in constraining the models, thus providing new information on 
both the geometry and physics of the HR 4049 circumstellar 
envelope detected by VINCI. This will indeed be possible with 
the advent of the MIDI and AMBER instruments and the aux- 
iliary telescopes to the VLTI facility in the next year or so. 
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Fig. 1. Best-fit through the observed visibility points of the normalized squared visibility curve provided by the simple models of 
HR 4049 described in the paper. The solid curve corresponds to the UD + Gaussian (star + envelope) model, the dashed-dotted 
line to the UD + UD model, both for a stellar to total flux ratio F=15%. The four panels refer to the different bins of azimuth 
(see Tabled}- 



